Spin- 1 ; Kondo effect in an InAs nanowire quantum dot: the Unitary limit, 
conductance scaling and Zeeman splitting 
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We report on a comprehensive study of spin-^ Kondo effect in a strongly coupled quantum dot 
realized in a high-quality InAs nanowire. The nanowire quantum dot is relatively symmetrically 
coupled to its two leads, so the Kondo effect reaches the unitary limit. The measured Kondo 
conductance demonstrates scaling with temperature, Zeeman magnetic field, and out-of-equilibrium 
bias. The suppression of the Kondo conductance with magnetic field is much stronger than would 
be expected based on a ^-factor extracted from Zeeman splitting of the Kondo peak. This may be 
related to strong spin-orbit coupling in InAs. 
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I. INTRODUCTION 



The Kondo effect^ is one of the most vivid manifesta- 
tions of many-body physics in condensed matter. First 
observed in 1930s in bulk metals through an anomalous 
increase in resistivity at low temperatures, it was later 
associated with the presence of a small amount of mag- 
netic impurities^ The modern theoretical understanding 
is that the single unpaired spin of the magnetic impu- 
rity forms a many-body state with conduction electrons 
of the host metal. This many-body state is characterized 
by a binding energy expressed as a Kondo temperature 
(Ik). When the temperature is decreased below Ik, the 
conduction electrons screen the magnetic impurity's un- 
paired spin, and the screening cloud increases the scatter- 
ing cross-section of the impurity. More recently, advances 
in microfabrication opened a new class of experimental 
objects semiconductor quantum dots in which a few elec- 
trons are localized between two closely spaced tunneling 
barriers 4 At the same time, it had been theoretically pre- 
dicted that an electron with unpaired spin localized in a 
quantum dot could be seen as an artificial magnetic im- 
purity and, in combination with the electrons of the leads, 
would display the Kondo effect 4»£ The first observation of 
Kondo effect in quantum dots was made in GaAs-based 
two-dimensional structures.— ~— Initially thought to be 
very difficult to observe in such experiments, the Kondo 
effect has now been seen in quantum dots based on a wide 
variety of nanomaterials such as carbon nanotubes , 11 ' 12 
Ceo molecules , 13 ' 14 organic molecules,— ~— and semicon- 
ductor nanowireSj^r— and has also been invoked to ex- 
plain behavior of quantum point contacts.— 

In this paper, we present a comprehensive study of 
the Kondo effect in a nanosystem of emerging interest, 
namely, InAs nanowircs grown by the vapor-liquid-solid 



(VLS) method.— Building on initial reports of Kondo 
effect in InAs nanowires ) 19 ' 20 we report Kondo valleys 
with conductance near 2e 2 /h in multiple devices and 
cooldowns. This high conductance, combined with tem- 
perature far below the Kondo temperature, allows quan- 
titative measurements of conductance scaling as a func- 
tion of temperature, bias, and magnetic field, which we 
compare to theoretical predictions independent of mate- 
rials system. The high (/-factor and small device area, 
characteristic of InAs nanowires, allows measurement of 
the splitting of the zero-bias anomaly over a broad range 
of magnetic field, and we find that splitting is pronounced 
at lower magnetic field than predicted theoretically. 



II. EXPERIMENT 

The quantum dot from which data are presented in 
this paper is based on a 50 nm-diameter InAs nanowire 
suspended over a predefined groove in a p + -Si/Si02 sub- 
strate and held in place by two Ni/Au (5nm/100nm) 
leads deposited on top of the nanowire. The leads' 450- 
nm separation defines the length of the quantum dot. 
The p + -Si substrate works as a backgate. The InAs 
nanowire was extracted from a forest of nanowires grown 
by molecular beam epitaxy on a (Oil) InAs substrate us- 
ing Au-catalyst droplets. Wires from this ensemble were 
found to have a pure wurtzite structure, with at most one 
stacking fault per wire, generally located within 1 fim 
from the tip. We therefore formed devices from sections 
of nanowire farther from the wires' end, with a reason- 
able presumption that the active area of each device is 
free of stacking faults. Schottky barriers, and screening 
of the electric field from the gate electrode by the source 
and drain electrodes, together create potential barriers 
next to the metal contacts. Thus electrons must tunnel 
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to the central part of the nanowirc (the quantum dot) 
and the contacts, giving rise to Coulomb blockade (CB). 
An SEM image of a typical device is shown in Fig. [lja). 
More details on growth, fabrication, and charging effects 
have been published previously^ 

Transport experiments were carried out in a dilution 
refrigerator with a base temperature Tb ase ~ 10 mK. 
All experimental wiring was heavily filtered and ther- 
mally anchored to achieve electron temperature close 
to cryostat base temperature, as verified in shot noise 
measurements^ Conductance measurements used stan- 
dard lock-in techniques with a home-built ultra-low-noise 
transimpedance preamplifier operated at frequencies of 
~ 2 kHz. Depending on the temperature T, the ac ex- 
citation bias was set in the range of 1-10 ^Vrms to keep 
it equal to or smaller than fc^T (ks is the Boltzmann 
constant). The magnetic field was applied perpendicu- 
lar to both the substrate and the axis of the nanowire. 
A schematic representation of the nanowire-based de- 
vice together with the experimental set-up is shown in 
Fig.HTb). 
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III. RESULTS AND DISCUSSION 

First we would like to outline the main features asso- 
ciated with the Kondo effect which were studied in our 
experiment. The conductance of a quantum dot weakly 
coupled to leads is dominated by CB, seen as nearly peri- 
odic peaks in the conductance as a function of gate volt- 
age, with the conductance strongly suppressed between 
peaks. Each peak signals a change in the dot occupancy 
by one electron. In contrast, a dot strongly coupled to 
the leads can show the Kondo effect, with the follow- 
ing signatures <£&2£ (1) the Kondo effect enhances con- 
ductance between alternate pairs of Coulomb blockade 
peaks (that is, for odd dot occupancy). These ranges of 
enhanced conductance are conventionally termed "Kondo 
valleys" . (2) Conductance in Kondo valleys is suppressed 
by increasing temperature. (3) Conductance in Kondo 
valleys is suppressed by applied source-drain bias (V s d), 
giving rise to a zero-bias anomaly (ZBA). The full width 
at half maximum (FWHM) of the zero-bias peak is of the 
order of 4fcsTk/e (e is the elementary charge). (4) In 
contrast to the conductance in the CB regime whose up- 
per limit is e 2 //i(2! the Kondo valley conductance can 
reach 2e 2 /h, equivalent to the conductance of a spin- 
degenerate ID wirei2£ In this limit, "valley" is a mis- 
nomer, as the valley is higher than the surrounding peaks! 
(5) The Kondo ZBA splits in magnetic field (B) with 
the distance between the peaks in bias being twice the 
Zeeman energy. (6) The dependence of the Kondo con- 
ductance on an external parameter A such as temper- 
ature, bias or magnetic field can be calculated in the 
low- and high-energy limits^ In the low-energy limit, 
I^bTk 3> A = {ksT, eV s d, |<?|/zbB}, the conductance has 



FIG. 1. (a) SEM image of a typical suspended nanowire- 
based quantum dot device used in the experiment. The scale 
bar corresponds to 1 (im. (b) Schematic representation of 
the nanowire-based quantum dot device and its experimen- 
tal setup, (c) The temperature dependence of the nanowire- 
based quantum dot conductance measured over a wide range 
of the backgate voltage V g . Five Kondo valleys are labeled 
I through V here. This identification of valleys will be used 
throughout the paper. Discontinuities in the temperature- 
dependence in Valley II are caused by device instability at 
this particular range of V g . (d) The gray-scale conductance 
plot in the V g — V s d plane measured in the same range of V 3 
as in (c) at temperature Tt aS e = 10 mK. Panels (a) and (b) 
are adapted with permission from A. V. Kretinin et at, Nano 
Lett. 10, 3439 (2010). Copyright © 2011 American Chemical 
Society. 



a characteristic quadratic Fermi-liquid behavior ; 14 i 30 — 



G(A) = Go 



(i) 



where Go = G(A = 0) and ca is a coefficient of order 
unity. Its numerical value is different for each parameter 
A, and depends on the definition of Tk- In the present 
paper, we use a convention^ used in many experimental 
papers and define Tk by the relation 



G(T = T K ) = 0.5G 



(2) 



In the opposite limit of high energy, when fceTk <C A, 
the conductance shows a logarithmic dependence. For 
example, as a function of tcmperature>i£ 



(3) 



G(T)ocG /ln z (^ 



There is no analytical expression for the intermediate 
regime, where the parameter A ss /c^Tk, but Numerical 



3 



Renormalization Group (NRG) calculations^ show that 
the connection between one limit and the other is smooth 
and monotonic, without any sharp feature at A = fc^TK- 

Before detailed consideration and discussion of the re- 
sults we give a broad overview of the experimental data 
used in this study. It will be followed by three subsections 
focusing on the observed unitary limit of the Kondo effect 
(Sec. IIII Aj) . conductance scaling with different external 
parameters (Sec. HiTB|) . and some peculiarities observed 
in the Zeeman splitting (Sec. IIII Cp . 

Figure [ljc) presents the linear conductance G as a 
function of the backgate voltage V g . Different color cor- 
responds to different temperature, ranging from 10 mK 
to 693 mK. The Kondo effect modifies the CB peaks 
so strongly that the separate peaks are no longer rec- 
ognizable and the simplest way to identify Kondo valleys 
is to look at the the gray-scale plot of differential con- 
ductance as a function of both V g and V s d ("diamond 
plot"), Fig. [Hd). Every Kondo valley is marked by a 
ZBA seen as a short horizontal line at V a d = 0. Different 
widths of ZBAs on the gray-scale plot reflect differences 
in the Kondo temperature. In these same Kondo val- 
leys, conductance decreases with increasing temperature 
(Fig. QJc)). Note that Kondo valleys alternate with val- 
leys having opposite temperature dependence or almost 
no temperature-dependence, corresponding to even oc- 
cupancy of the quantum dot. A small unnumbered peak 
at about V g = -2.95 V departs from the general pattern 
of conductance observed in the experiment. Most likely, 
this feature, which occurs for even occupancy, is associ- 
ated with transition to a triplet ground state, and thus 
emergence of spin- 1 and singlet-triplet Kondo effect M-~— 
However, it is difficult to conclusively identify the nature 
of this anomaly since its temperature and bias dependen- 
cies are weak. 

All conductance peaks shown in Fig.QJc) exceed e 2 /h, 
reflecting Kondo-enhanced conductance and relatively 
symmetric coupling to the two leads. In particular, con- 
ductance around V g = -3.1 V in Valley III reaches the 
unitary limit of 2e 2 /h, to within our experimental accu- 
racy. 



A. Kondo effect in the unitary limit 

To realize maximum conductance in resonant tunnel- 
ing, the quantum dot should be symmetrically coupled 
to the leads. In the conventional case of CB, electro- 
static charging allows only one spin at a time to tunnel, 
limiting the maximum conductance through the dot to 
e 2 //ii2! The Kondo effect dramatically changes the sit- 
uation by forming a spin-degenerate many-body singlet 
state, enabling both spins to participate in transport in 
parallel so that Kondo conductance can reach its unitary 
limit at 2e 2 /h££ Experimentally, the unitary limit, first 
observed by van der Wiel et al££- in a GaAs-based gate- 
defined quantum dot, remains the exception rather than 
the rule, because it requires being far below the Kondo 




-3.1 -3.05 
V V 

9 V 1 



FIG. 2. The Kondo effect in its unitary limit. The main plot 
shows the linear conductance G in Valley III, as a function of 
backgate voltage V g at different temperatures. The dark blue 
curve corresponds to the lowest temperature of 10 mK. Inset: 
The red triangles correspond to the temperature dependence 
of the conductance at a fixed V g = -3.107 V (marked by the 
red triangle in the main graph). The blue curve represents 
the result of approximation with Eq.Q where Go = 1.98e 2 //t 
and T K = 1.65 K. 



temperature, having symmetric tunnel coupling to the 
two leads, and having precisely integer dot occupancy. 

Figure [2] presents a zoomed-in view of Valley III 
from Fig. [IJc), showing the Kondo effect in the unitary 
limit. Note how the conductance maximum gradually ap- 
proaches 2e 2 /h with decreasing temperature. Here the 
limit is reached only at some particular V g , showing a 
peak instead of an extended plateau as reported by van 
der Wiel et al.££- Since tunneling is so strong that level 
widths are almost as large as the Coulomb interaction 
on the dot, the dot occupancy rid is not well-quantized 
but rather changes monotonically, passing through rid = 
1 (n-f = rii = V2) a t Vg ~ —3.1V, where the unitary 
limit is observed. In accordance with the Friedel sum 
rule the conductance of the dot is predicted to depend on 
the dot occupancy nf,j. as G(t, 4-) = (e 2 /^-) sin (7rnf ,,[.). 
So the sum of the conductances is 2e 2 /h when rid = 1- 
Note that the Kondo conductance shown in Fig. [IJc) al- 
ways exceeds 1.3 e 2 /h for different dot occupancies, show- 
ing that the wave function overlap with the two leads is 
rather equal: the two couplings are within a factor of 
four of each other over this whole range, suggesting that 
disorder along the nanowire and especially at the tunnel 
barriers is quite weak. To extract the Kondo temperature 
we apply a widely used phcnomcnological expression^ for 
the conductance G as a function of temperature: 



G(T) = G 1 + {T/T' K ) 



(4) 



where Go is the zero-temperature conductance, T' K = 
T K /(2 1 / S - l) 1 / 2 , and the parameter s = 0.22 was found 
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to give the best approximation to NRG calculations for 
a spin- A Kondo system^ Here, the definition of Tk is 
such that G(T K ) = G /2. The inset of Fig. [2] shows the 
conductance for different temperatures at V g = -3.107 V 
(marked by the red triangle in the main figure). The blue 
curve in the inset represents the result of the data approx- 
imation using Eq. (j4|) where the fitting parameters Go and 
T K are (1.98±0.02)e 2 /^ and 1.65±0.03 K^l respectively, 
showing that the system is in the "zero-temperature" 
limit at base temperature, T^/T^ asc 165. 



B. Conductance scaling with temperature, 
magnetic field and bias 

As noted above, the Kondo conductance as a function 
of temperature, bias or magnetic field should be describ- 
able by three universal functions common for any sys- 
tem exhibiting the Kondo effect. Before discussing ex- 
pectations for universal scaling we describe in detail how 
temperature, magnetic field and bias affect the Kondo 
conductance in our experimental system. 



1. Kondo conductance and Kondo temperature at zero 
magnetic field 

For a more detailed look at the spin- A Kondo effect 
at B = 0, we select the two Kondo valleys IV and V 
(see Fig. He)). The zoomed-in plot of these two val- 
leys is shown in Fig. |2Ka,b). The coupling to the leads, 
and hence the Kondo temperature, is much larger in Val- 
ley V than in Valley IV. Valley IV shows a typical ex- 
ample of how two wide Coulomb blockade peaks merge 
into one Kondo valley as the temperature decreases be- 
low Tk^&2£ Valley V, in contrast, does not evolve into 
separate CB peaks even at our highest measurement tem- 
perature of 620 mK. Also, as seen from Fig. EKb), the 
width of the ZBA, which is proportional to Tk, is larger 
for Valley V. To illustrate this, in Fig. |U[a,b), we plot 
the conductance as a function of V sc i at different temper- 
atures for two values of V g (marked by red triangles in 
Fig. EKa)) corresponding to the two valleys. In addition 
to the ZBA of Valley IV being significantly narrower than 
that of Valley V, at the highest temperatures, the ZBA of 
Valley IV is completely absent, while the ZBA of Valley 
V is still visible, pointing to a significant difference in Tk- 
To quantify this observation we found Tk as a function of 
V g for both valleys by fitting the temperature-dependent 
conductance using Eq. The result of this fit is pre- 
sented in Fig. |U(c,d). Tk shows a parabolic evolution 
across each Valley, with Tk ranging from 0.3 K to 1 K 
for Valley IV and from 1.3 K to 3 K for Valley V. This 
significant difference in Tk correlates with the difference 
in the ZBA width shown in Fig. Hta,b). However, the 
relation between the FWHM of the ZBA peak and Tk is 
more ambiguous due to out-of-equilibrium physics^ 

To understand the dependence of Tk on V g and to 
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FIG. 3. (a) The detailed measurement of the conductance 
temperature dependence shown in Fig. QJc), Valleys IV and 
V. The red triangles mark two values of V g = -2.835 V and 



V a 



-2.680 V for which the conductance as a function of 



Vsd is plotted in Fig. |4ja) and Fig. [U[b), respectively, (b) 
The gray-scale conductance plot in the V g — V 3 d plane was 
measured in the same range of V g as in (a), at temperature 
T = 10 mK. 



extract some relevant parameters of the system we use 
an analytic prediction for the dependence of the Kondo 
temperature based on the microscopic parameters in the 
Kondo regime of the single-impurity Anderson models 



'TU 



exp 



7reo(eo+{7) 
TU 



(5) 



Here T is the width of the resonant tunneling peak, 
U = e 2 jGtot is the charging energy (C tot is the total 
capacitance of the dot), and Eq is the energy of the res- 
onant level relative to the Fermi level. As Tk is derived 
from the conductance [c.f. text following Eq. (0J)], the 
prefactor J? NRG in Eq. ([5]) of order unity was calibrated 
using the NRG. To this end, we calculated the conduc- 
tance G(T) for the single-impurity Anderson model at 
eo = —U/2, for fixed U/T ~ 4.5. The requirement that 
G(T = T K )/G(0) = G /2 fixes the prefactor in T K to 
7y NRG ~ 1.10, which we took constant throughout. ry NRG 
does vary as a function of U /T within a few tens of per- 
cent, due to the exponential sensitivity of Eq. ([5]), how- 
ever, since U and T are already pretty well-constrained in 
our case, this results in negligible variations in our fitted 
U, e or T. 

To determine the parameters U, £q and T, we proceed 
as follows. The value of U ~ 400/ieV was found from 
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FIG. 4. Nonlinear conductance as a function of V s d around 
zero bias for different temperatures at V a = -2.835 V (a) and 
V g — -2.680 V (b), near the centers of Kondo valleys IV and 
V. The color scale is as in Fig. |3ja). (c,d) The Kondo tem- 
perature 7k, plotted on a semi- log scale, as a function of V g 
for these same valleys. Panel (c) corresponds to Valley IV 
and panel (d) to Valley V. Blue curves in both panels show 
fits of Eq. © to data, with T IV ~ Y7Q^eV for Valley IV and 
Y v « 435/xeV for Valley V. 



Fig. 02[b) for Valley IV (we assume the value is equal 
for Valley V, though it may be slightly lower, given the 
stronger tunnel coupling there). To relate £o and V g , 
we used a simple linear relation V g — Vgo = «£o with 
the lever arm a = Ctot/G g , where V^o is the position 
of the Coulomb peak and C g is the gate capacitance. 
Here C to t = e 2 /U and C g = e/AV g where AV^ is the 
CB period. T was determined by fitting the curvature of 
log Tk with respect to gate voltage in Fig.|U[c,d), yielding 
T JV k 176 neV and T v re 435 ^eV for Valley IV and V, 
respectively. 

As noted above, the predicted dependence of Tk in 
Eq. ((SJ is based on the Anderson model in the Kondo 
regime (e /T < -1/2)M The fitting of the data with 
Eq. (|5|), however, gave Eo/Tiy ~ —1.1 and So/Ty ~ —0.5 
in the centers of Valley IV and V, respectively. So the 
Kondo regime {|eo|, |eo + U\} > Y/2 is reached only near 
the center of Valley IV and only at the very center of 
Valley V. The rest of the gate voltage range in these 
Valleys is the mixed valence regime, where charge fluc- 
tuations are important and Kondo scaling should not be 
quantitatively accurate^ Note: our NRG calculations 
show that the deviations from universal scaling up to 
£o ~ — r/2 should be small for T < Tk- In any case, we 
have not attempted to take into account multiple levels 
in our calculations, which could quantitatively but not 



qualitatively modify the predicted behaviors. 



2. Kondo conductance at non-zero magnetic field 

The Kondo effect in quantum dots at non-zero mag- 
netic field is predicted and observed to exhibit a Zeeman 
splitting of the ZBA by an energy A = 2|g|^B-B&£ (g is 
the g- factor, /xb is the Bohr magneton), which is a direct 
consequence of the (now-broken) spin-degeneracy of the 
many-body Kondo singleti 41 ' 42 

To analyze the Zeeman splitting in our nanowire-based 
quantum dot we focus on Kondo Valley IV. The Kondo 
ZBA at zero field, seen in a zoom- in in Fig. [SJa), is sup- 
pressed at B = 100 mT, but recovers once a bias of 
~ 40 /j,V is applied (Fig.[SJb)). Contrary to earlier obser- 
vations in InAs nanowires^ we find that the g-factor at a 
given field is independent of V g as illustrated by the par- 
allel slit-like shape of the Zeeman splitting (Fig. [5jb))4£ 
The gray-scale conductance plot in Fig. [5f c) presents the 
evolution of the Zeeman splitting with magnetic field at 
fixed V g = -2.835 V, marked by the cross in Fig.[5ja) (for 
the associated ZBA measured at B = refer to Fig.@|a)). 
The plot shows the splitting in bias A/e to be almost 
linear in magnetic field, which allows us to deduce the 
value of the (/-factor by fitting the data with a linear de- 
pendence V sd = ±\g\fi B B/e for 30 mT< B < 100 mT. 
Two red lines in Fig. [5tc) show the result of fitting with 
\g\ = 7.5±0.2 (the meaning of the dotted green lines will 
be discussed below). This number is smaller by a factor 
of two than the InAs bulk value of \g\ = 15, possibly due 
to the reduced dimensionality of the nanowire device* 4 ^ 
and it is consistent with previous measurements J£ 

We now compare the dependence of the Kondo con- 
ductance on the temperature and magnetic field, respec- 
tively. In order to do so we plot on the same graph 
G{T,B = 0) and G(T = T hasc ,B) both taken in equi- 
librium at V g = -2.835 V (Fig. [Sfd)). In order to quan- 
titatively compare the effect of magnetic field to that of 
temperature, we associate each magnetic field value with 
an effective temperature Tb(B) = \g\^B / k B . where 
\g\ = 7.5 is extracted from the linear Zeeman splitting of 
peaks in differential conductance. The comparison of the 
linear conductance data is presented in Fig. [5^d) , where 
G(T) is shown by the blue squares, G{B) by the red 
triangles. In this same plot, theoretical predictions are 
shown as curves: blue for G(T) and red for G(B). Note 
that for |<7| = 7.5 (this value extracted from the split- 
ting of the differential conductance peaks) , the blue and 
red curves differ substantially for essentially all nonzero 
values of their arguments, with magnetic field having a 
much weaker predicted effect than temperature. Sur- 
prisingly, in light of this theoretical prediction, the two 
sets of experimental data lie almost on top of one an- 
other up to about 200 mK « T K . The NRG results for 
G(T = 0, i?) 42 ' 45 have been checked against exact Bcthe 
ansatz calculations 4 ^ 4 ^ for G(T = 0, B) (dashed black 
curve in Fig. [Sfd)) and are seen to be in excellent agree- 
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of magnetic field effects on linear and differential conduc- 
tance is a conundrum. We hope this work will stimulate 
further theory and experiment to address this issue. 



3. Universal conductance scaling 
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FIG. 5. The Zeeman splitting of the Kondo ZBA measured 
at T — 10 mK. (a) The gray-scale conductance plot of Kondo 
valley IV (see Fig.[3ja)) measured at B — 0. (b) The same as 
in (a) but at B = 100 mT. (c) Gray-scale conductance plot in 
the V s d — B plane measured at fixed V g = -2.835 V denoted 
by the cross in panel (a). The red dashed lines represent the 
result of the fitting with expression V s d = ±\g\fJ,BB/e, where 
\g\ = 7.5±0.2. Vertical blue dashed line marks magnetic field 
value 0.5fcsTK/|(/|^iB as a reference for the onset of Zeeman 
splitting (here Tk = 300 mK). While \g\ = 7.5 gives the best 
match to linear Zeeman splitting, \g\ = 18 (green dotted lines) 
could account for the fact that Zeeman splitting is resolved at 
very low field, (d) Conductance at V s d = as a function of T 
(blue squares) and as a function of the effective temperature 
Tb = \g\nsB / ks (red triangles). The solid blue curve shows 
G{T) from NRG, the solid red curve G(B) from NRG, and 
the dashed black curve G(B) from exact Bethe Ansatz (BA) 
calculations for the Kondo mode h 45 ' 46 These assume \g\ — 7.5. 
For NRG and BA calculations of magnetic field dependence, 
additional curves (solid green and dashed brown) are plotted 
for \g\ = 18, showing better match to linear conductance data 
- though not to the differential conductance in (c) above. 



ment, so the disagreement between theory and experi- 
ment is not related to a particular calculational frame- 
work. Were we to assume \g\ — 18, we could explain the 
experimental magnetic field dependence of linear conduc- 
tance G(T — 0,B), as shown by alternative curves (solid 
green and dashed brown) plotted in [SJd). This value of 
g is within the realm of possibility for InAs nanowires. 20 
However, we are inclined to rely on the g value of 7.5 ex- 
tracted from the splitting of the peaks in the differential 
conductance. With \g\ — 18 we would have the puzzling 
result that the splitting of peaks in differential conduc- 
tance would be less than half the expected 2\g\fisB (see 
dotted green lines in Fig. 03c)), which would be hard to 
explain. Regardless, the mismatch between the strength 



In testing universal conductance scaling, we concen- 
trate first on the scaling of the linear conductance with 
T and B. In the case of temperature dependence, the 
universal scaling function has the form of Eq. Q. This 
expression has been applied to a wide variety of exper- 
imental Kondo systems^ ' 11 ! 14 ' 19 and after expansion in 
the low-energy limit (T/Tk <C 1) it becomes Eq. ^ de- 
scribing the quadratic dependence on temperature^ 2 - 



G W Go 



1 - c T (T/T K )' 



(6) 



where c T = c A = s(2 1/s - 1 ) = 4.92 and s = 0.22 is taken 
from Eq. ([4]). Note that this coefficient ct is about 10% 
smaller than the more reliable value ct ~ 5.3 S 30 ' 33 ' 47 ' 48 
found from the NRG calculations on which the phe- 
nomenological form of Eq. (Ql is based (see Table HJ 42 
Since Eq. ([4]) is independent of the particular sys- 
tem, it can be used as the universal scaling function 
G/Gq = /(T/Tk). Figures [HJa,b) show the equilibrium 
Kondo conductance (1 — G/Gq) of valleys IV and V (see 
Fig. [3^a)) plotted as a function of T/Tk, taken at dif- 
ferent V g . Here, the values of Go and Tk are found by 
fitting the data with Eq. © for T < 200mK (for higher 
temperatures the conductance starts to deviate from the 
expected dependence due to additional high-temperature 
transport mechanisms). As seen in Fig. Oa,b), all the 
data collapse onto the same theoretical curve (dashed) 
regardless of the values of V g or Tk- In the low-energy 
limit T/Tk < 0.1 the conductance follows a quadratic de- 
pendence set by Eq. (fTJ) with coefficient ca = ct = 5.38 
as shown by the dotted line. As noted above, in the low- 
energy limit, the phenomenological expression Eq. (|4]) 
is less accurate and shows a quadratic dependence with 
ct = 4.92. This explains why the dashed and dotted 
curves in Fig.[S)Ia,b) do not coincide at T/Tk < 0.1. 

It should also be possible to scale G(B) as a function of 
a single parameter Tb/Tk. As an example, we present in 
Fig.EHa) scaled G(B) data from Fig.^d). At low fields, 
the measured conductance is found to depend on B ac- 
cording to Eq. ([T]), with the coefficient ca = cb ~ &r- 
This equality has also been independently checked by fit- 
ting the G(B) and G(T) data for T/T Kl T B /T K < 0.1 
with Eq. (HQ). The ratio between the two fit coefficients, 
Cb/ct, is approximately 1 (cb/ct = 0.95 ± 0.2), strongly 
counter to the theoretical expectations where cb ~ 0.55 48 
and cb/ct ~ 0.101 (see Table [J. To illustrate this dis- 
crepancy, we plot Eq. (JT|) with ca~cb= 0.55 in Fig^a) 
(dash-dot). The reason for such a dramatic difference 
in G(B) dependence between theory and experiment for 
both low- and intermediate-field range is unclear. We 
speculate that the spin-orbit interaction, previously ob- 
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FIG. 6. (a,b) The equilibrium conductance of Kondo valleys 
IV (a) and V (b) at different V g , scaled as a function of a sin- 
gle argument T/Tk (blue squares) and Tb /Tk (red triangles), 
where Tb = \g\^B /k^.. The dashed curve shows the univer- 
sal function described by Eq. Q. The dotted line represents 
the low-energy limit of Eq. ([T]) with ca = ct = 5.38. The 
dash-dotted line shows the theoretically predicted low-field 
scaling of G(B) with cb = 0.55. The values of Go and Tk were 
found by fitting the data with Eq. (0}, see Sec. IIIIB II For 
values of V g refer to Fig.[4jc,d) and Fig.[5jd). (c,d) The scaled 
conductance AG/5 = (1 - G(T, V sd )/G(T,0))/a, where & = 
CTa/(l + c T (7/a-l))(T/Tk) 2 , versus (eV s d/k B T K ) 2 taken at 
several V g along Kondo valleys IV (c) and V (d). For valley IV 
the backgate voltage was chosen from the range V g — -2.82 V 
to -2.85 V with 5 mV step and for valley V from the range 
V g = -2.68 V to -2.72 V with 20 mV step. Different colors 
of the data points represent different temperatures (9.5 mK, 
12.9 mK, 22.4 mK, 32.6 mK, 46.1 mK, 54.2 mK). The dashed 
line shows the corresponding scaling function given by Eq. © 
with a = 0.18 and 7 = 1.65. 



served in InAs nanowire-based quantum dots^, may play 
a role. 

It is important to note that in order for the universal 
scaling G(B) to be valid, the coefficient Go in Eqs. ([1]) 
and J3|) should be independent of B. In the case of GaAs 
quantum dotsL ^ 26 ' 51 with \gGaAs\ = 0.44 the magnetic 
field required to resolve the Zeeman splitting is high and 
the orbital effects of that field contribute significantly, 
resulting in a B-dcpcndent Go, even for a field paral- 
lel to the plane of the hetcrostructure. In contrast, in 
our InAs nanowire-based quantum dot, with large g- 
factor and small dot area S = 50 nm x 450 ran, Kondo 
resonances are suppressed (split to finite bias) at fields 
smaller than that required to thread one magnetic flux 
quantum B < {h/e)/S ~ 180 mT, thus making the or- 
bital effects negligible and Go magnetic field indepen- 



dent. 

Now that the scaling of the linear conductance has 
been established, including the stronger-than-expected 
effect of magnetic field, we examine how the out-of- 
equilibrium conductance scales as a function of bias and 
temperature G/G = f(T/T K ,eV sd /k B T K ). The func- 
tion used to test the universal scaling in a GaAs quan- 
tum dot^ and in a single- molecule device^ originates 
from the low-bias expansion of the Kondo local density 
of states^ and has the following form: 



G(T,V sd )=G(T,0) 



eVsd 
k B T K 



(7) 

The coefficients a and 7 relate to the zero-temperature 
width and the temperature-broadening of the Kondo 
ZBA, respectively. The zero-bias conductance G(T, 0) 
is defined by Eq. (|6|). The coefficients a and 7 are in- 
dependent of the definition of the Kondo temperature 
and in the low-energy limit Eq. ([7]) reduces to the theo- 
retically predicted expression for non-equilibrium Kondo 
conductance^ 



G(T,V sd )-G(T,0) 
c T G 



«( 



eV sd 



ct7 



(*)'( 



(8) 

The independence of a and 7 on the definition of 
Kondo temperature is important: though we have chosen 
an explicit definition for Tk, consistent with the choice 
used for most quantum dot experiments and NRG calcu- 
lations, other definitions may differ by a constant multi- 
plicative factor. 

Figures |U[c,d) show the scaled finite-bias conductance 
(1-G(T, V sd )/G(T,0))/a, where 5 = c T a/{l + c T {i/a- 
I))(T/Tk) 2 , versus (eKd/^sTx) 2 , measured at different 
temperatures and a few values of V g . The conductance 
data are fit with Eq. [7] using a procedure described by 
M. Grobis et alM> with two fitting parameters a and 7. 
The range of temperatures and biases used for the fit- 
ting procedure was chosen to be close to the low-energy 
limit, namely T/T K < 0.2 and eV sd /k B T K <0.2, which 
is comparable to the ranges used in Ref. l32t Averag- 
ing over different points in V g gives a = 0.I8±0.0I5 and 
7 = 1 .65±0.2 for Valley IV. Despite Valley V being in the 
mixed- valence regime, the parameters a and 7 are close 
to those found for Valley IV. The scaled conductance in 
both cases collapses onto the same curve, shown by the 
dashed line, for ± {eV s d / k^T^) 2 < 0.1, though the data 
from Valley V deviate more from the predicted scaling. 
This is not surprising because the Valley V data are in 
the mixed-valence regime, beside that the bias can cause 
additional conduction mechanisms due to proximity of 
the Coulomb blockade peaks. 

Overall, the value of a obtained in our experiment 
is larger than previously observed in a GaAs do t 32 ' 53 
(a = 0.1) and single molecule^ (a = 0.05). The exact 
reason for this discrepancy is unknown, but the smaller 
ratio Tbasc/Tk may play a role. 
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There is a large number of theoretical works de- 
voted to the universal behavior of finite-bias Kondo 
conductance based on both the Anderso n 33 ! 48 ' 54 " — and 
Kond o 29 ' 31 ' 59 " — models. Early predictions based on an 
exactly solvable point of the anisotropic nonequilibrium 
Kondo mode l 31 ' 59 ' 60 yielded a value a = cv/ct = 3/tt 2 ~ 
0.304. This turned out to be in disagreement with ex- 
periment, which is not surprising, since this coefficient 
is not universal and hence will not be the same for the 
isotropic Kondo models. A number of subsequent papers 
that used a Fermi-liquid approach to treat the strong- 
coupling fixed point of the Kondo model 29 ' 63 or studied 
the U — s- oo limit of the symmetric Anderson modelj^ 4 -"— 
all found a = 3/(2ir 2 ) sa 0.152. Our measured value of 
a = 0.18 is in a good agreement with this prediction. A 
Bethe-Ansatz treatment of the nonequilibrium Anderson 
model^ 8 - yielded a different result, a — 4/7r 2 , but this was 
obtained using some approximations and was not claimed 
to be exact. 

Some of the more recent theoretical papers have stud- 
ied the ay-coefficients for the nonequilibrium Anderson 
model under less restrictive conditions, i.e., allow for a 
left-right asymmetry and a non-infinite U , in an attempt 
to explain the experimental results of Refs. [l4| and [32|. 
J. Rincon and co-authors^ 5 - found that by setting U to 
be finite the expected value of a is decreased from 0.152 
to 0.1, but 7 remains ~ 0.5. Later, P. Roura-Bas 5 ^ came 
to a similar conclusion considering the Anderson model 
in the strong coupling limit in both the Kondo and the 
mixed-valence regimes. It was shown^ 6 - that a reduces 
from 0.16 to 0.11 if some charge fluctuation is allowed by 
shifting from the Kondo to the mixed- valence regime, and 
the parameter 7 is not necessarily temperature indepen- 
dent. In an attempt to explain the small a observed in 
molecular devices 1 ^ Sela and MaleckiS evaluated a model 
for the Anderson impurity asymmetrically coupled to the 
leads. They concluded that deep in the Kondo regime a 
takes the value of 3/(27r 2 ) m 0.152 independent of cou- 
pling asymmetry. However, if U is made finite or, in 
other words, some charge fluctuations are included, the 
parameter can vary within the range 3/(47r 2 ) < a < 3/tt 2 
(0.075 < a < 0.3) depending on the asymmetry of the 
tunneling barriers. Despite the fact that our system is far 
from the strong coupling limit (U ~ T, instead of U 3> T, 
see Sec. IIII B~T|) . the observed value of a = 0.18 is a good 
match to the strong-coupling prediction. 

From temperature, magnetic field, and bias scaling of 
the measured conductance, we are able to define a com- 
plete set of coefficients ca to be used in Eq. ([TJ in order 
to describe the Kondo effect in the low-energy limit: 

G(T) = G [l-c T (T/T K ) 2 ], 

G[B) = G Q [1 - c B (\g\fx B B/k B T K ) 2 ], 

G{V sd ) = G [l - c v (eV sd /k B T K )% 

where Go is the conductance at zero temperature, mag- 
netic field and bias, ct « 5.6±1.2, cb ~ 5. 1±1 .1, 



cy = ctol w 1.01±0.27. The substantial uncertainties 
originate from the small number of experimental points 
satisfying the requirement of low temperature, field, and 
bias used during fitting with Eq. ([1]). Table U summa- 
rizes the experimental value of these three parameters 
and compares to their theoretical predictions. (The pa- 
rameter a discussed above is denoted by ay in the table.) 
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FIG. 7. (a) The non-equilibrium Kondo conductance as a 
function of V s d for several values of B (open blue squares). 
The solid red curves represent the approximation of the 
data made with the sum of two Fano-shaped peaks and 
a cubic background. (b) The normalized Zeeman split- 
ting A/[2|<7|/^b-B] as a function of B data acquired from 
the peak maximum search (blue squares) and after fitting 
with two asymmetric peak shapes (red triangles). The 
vertical blue and green dashed lines denote magnetic field 
of 0.5fcBTk/|(?|A 1 B and fcsTk/lslA'B correspondingly (here 
\g\ = 7.5 and T K = 300 mK). 



C. Zeeman splitting 

At non-zero magnetic field, the spin degeneracy of 
the Kondo singlet is lifted and the linear conductance 
through the dot is suppressed.— To recover strong trans- 
port through the dot a bias of i^A/e = ±|<?|/isi?/e 
should be applied in order to compensate for the spin- 
flip energy. As a result, in experiments, the ZBA is 
split into two peaks separated by eA = 2\g\fiBB/efi& 
providing information on the effective g-factor. This is 
why the splitting of the Kondo conductance feature has 
become a popular tool for evaluating the value and be- 
havior of the g-factor in quantum dots made of differ- 
ent materials ! 12 ' 16 ' 17 ' 19 ' 20 ' 26 ' 68 In this section, we discuss 
two unexpected features related to the Zeeman splitting. 
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TABLE I. Summary of theoretically predicted parameters ct, cy, cb, and B c and their experimental values. The second 
column lists the values of the parameters c' A appearing in G(A) = Go[l — c^(A/fcBTb) 2 ], using a definition for the Kondo 
scale that is widespread in theoretical papers, namely, To — l/(4xo), where xo is the static impurity spin susceptibility at 
T — 0. This definition of the Kondo temperature differs from the Tk used in this paper, i.e., G(Tk) = G(0)/2, by the factor 
Tk/To = 0.94. 64 Thus, the coefficients c A defined in our Eq. (1) and listed in the fourth column are related to those in the 
second by ca/c' a = (Tk/To) 2 . We cite only references that are relevant for the symmetric Anderson model in the large- U limit, 
where the local occupancy is 1; generalizations for the asymmetric Anderson model may be found in Refs. |55|. l57l 15a. and 163. 
The last row lists values for the critical magnetic field B c beyond which the Kondo ZBA splits and it is expressed in units of 
Tk defined by Eq. ^ (Theory: column 2; Experiment: column 5). 



Parameter 


Predicted c' A 


Ola = 


Ca/c t 


ca = c' A (Tk/To) 2 


Experimental value 


ct 

CV 

cb 


7T 4 /16 « 6.08822,30,33,47,48,63,65 

3tt 2 /32 a 0.9253Ufr&iS3iS2 
7T 2 /16 « 0.617 29 ^^ 


3/(2rr 2 ; 
1/rr 2 : 


1 

) w 0.152 
a 0.101 


5.38 
0.82 
0.55 


5.6±1.2^ 

1.01±0.27 a ; 0.670^^ 0.304^ 
5.1±l.l a 


IglfJ-sBc/ksTK 


i.06fl;i.04aa ; i.itt 








<0.5 a ; 0.5^; l^i; 1.5*2 



a Present experiment 



First, the minimal value of field needed to resolve the 
Zeeman splitting is lower than expected. Second, the 
splitting is weakly sublinear with magnetic field at larger 
fields. 

Some attention has been previously paid to the value 
of the critical field B c at which the splitting of the 
Kondo ZBA occurs. The theory developed by one of 
the present authors^ 2 - predicts the value of the criti- 
cal field at T/T K < 0.25 to bc B c = lMk B T K /\g\^ B , 
with similar values being found by other authors i 66 ' 67 ' 69 
Treating nonequilibrium more realistically gives a slightly 
larger value.— Recent work by the authors, using den- 
sity matrix approaches ) 70 ' 71 suggests that a precise de- 
termination of the critical field is a numerically diffi- 
cult task, which will require further work in order to 
establish this beyond any doubt. There are also some- 
what conflicting experimental data on this issue. The 
value of B c predicted by Costi^ and Hewson et alf& 
seems to agree with the experimental findings for GaAs 
dotsj^i however, in gold break junctions^ the onset of 
the splitting was measured at O-SfceTk/I^I^B and in the 
case of carbon nanotubes±£ at about l.bk^T^/\g\^. In 
our case, Tk = 300 mK (see Fig. Hfc)), thus the pre- 
dicted B£2- is expected to be ~ 60 mT (for \g\ = 7.5), 
more than twice as large as that observed experimen- 
tally: as seen in Fig. [TJa) and Fig. [Sfc) , the splitting is 
already well-resolved at B = 30 mT, which corresponds 
to 0.5/cBXk/|g|/XB, ^ ne same as the result for gold break 
junctions^ Such a wide deviation of B c found for vari- 
ous Kondo systems (see Table HJ may be associated with 
a different width of ZBA (relative to Tk) in the vari- 
ous experiments. Since the conductance peak discussed 
here (see Fig. 0Ja)) is rather narrow, most likely due to 
the relatively low temperature T/Tk_ ~ 1/30, it is possi- 
ble to resolve the splitting onset at lower magnetic field. 
The analysis of the non-equilibrium scaling parameters, 
described in Sec. IIIIB 3[ confirms the above assumption. 

Finally, we discuss the evolution of the splitting A 



with magnetic field. Theory predicts that the peaks in 
the spectral function for spin-up and spin-down electrons 
should cling closer to zero energy at relatively low mag- 
netic fields than might naively be expected, so that A 
should bc suppressed by up to « 1/3 in the low-field 
limiti 47 ' 72 " — One recent experimental report corrobo- 
rates this predicted trend of suppressed splitting at low 
field.— But the variety of deviations from linear splitting 
in experiments - especially near the onset of splitting - 
is large i 12 ' 51 To make small variations in A more visible, 
we plotted the normalized value 8(B) = A/[2|<?|/j,b-B] in 
Fig. [T^b). The value of A was deduced from a simple 
peak maximum search (blue squares) and by fitting the 
data with the sum of two asymmetric peak shapes and 
some background (red triangles). To fit G as a function 
of V s d we used a combination of two Fano-shape asym- 
metric peaks on a cubic background: 

G(V sd ) = A, Lp=g^ +A 2 [^Ij +B\V sd \* + C. 

Here A\ and A 2 are the amplitudes, T\ and T 2 are the 
widths, gi and q 2 are the asymmetry parameters of the 
two Fano resonances positioned at dc bias Vi and V2, 
respectively. Parameters B and C characterize the cu- 
bic conductance background. Without the cubic back- 
ground, the positions of the conductance peaks, which 
correspond to Fano resonances at V\ and V 2 would be 
V p i = Vi + Ti/qi and V p2 = V 2 + T 2 /q 2 . The peak 
separation is deduced from the fit according to the equa- 
tion A/e = V p2 — V p i- The quality of this fit is shown 
in Fig. Oa) by red solid curves. It is clear that at 
B > 100 mT the splitting is sublinear in magnetic field. 
Coincidence of the splitting data extracted by two differ- 
ent methods (blue triangles and red squares in Fig.|7£b)) 
makes us believe that this effect is genuine and not an 
artifact due to weakly bias-dependent background con- 
ductance. In contrast, splitting extracted from our data 
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at low fields B < fcBlk/|<7|A t B is dependent on the extrac- 
tion method used, so we do not wish to make quantitative 
claims for the magnitude of splitting in that field range. 
Our results differ from previous observations mainly in 
that a sublinear field splitting occurs also at higher fields 
and not only at the onset of the splitting . 12 ' 51 We are 
unaware of any theoretical predictions which would ex- 
plain such sublinear splitting or effective reduction in the 
g-factor at higher fields. 

Previous theoretical works on the Kondo model pre- 
dicted a suppressed splitting S(B) = A/2\g\[i B B in- 
creasing monotonically toward one for gfisB 3> fes^k 
with logarithmic corrections , 72 1 76 1 77 For the Anderson 
model, similar results have been found with 5(B) ris- 
ing monotonically with increasing 5 . 73 ' 78 ' 79 However, in 
some works ; 67 ' 73 ' 78 5(B ^> UbTk) is found to exceed one, 
whereas in other work s 47 ' 79 S(B 3> fc^Tx) remains be- 
low one. This discrepancy between different approaches 
is likely due to different approximations and the extent 
to which universal aspects as opposed to non-universal 
aspects are being addressed and remains to be clari- 
fied. For example, it is known that extracting peak po- 
sitions in equilibrium spectral functions within NRG is 
problcmatici 67 ' 79 ' 80 Extracting a Zeeman splitting from 
experimental dI/dV s d at finite bias and large magnetic 
fields is also complicated by the increasing importance of 
higher levels and non-equilibrium charge fluctuations^ 
Nevertheless, our results for 8(B » fcfiTk) in Fig. [7jb) 
exhibit a monotonically decreasing 8(B) in the high field 
limit for B > 1.5fcBXk/|<?|MB- This contrasts to current 
theoretical predictions. As we cannot exclude the contri- 
bution of orbital effects at higher B, the magnetic fields 
used to determine the (/-factor were chosen to be smaller 
than 100 mT (flux through dot < 0.6$ ). 

IV. CONCLUSION 

In conclusion we have performed a comprehensive 
study of the spin-i Kondo effect in an InAs nanowire- 
base quantum dot. This experimental realization of a 
quantum dot allowed us to observe and thoroughly ex- 
amine the main features of the Kondo effect including 
the unitary limit of conductance and dependence of the 
Kondo temperature on the parameters of the quantum 
dot. Also the Kondo temperature's quantitative relation 
to the Kondo ZBA shape, Zeeman splitting of the ZBA 
and scaling rules for equilibrium and non-equilibrium 
Kondo transport were studied. A previously undetected 
dependence of the (/-factor on magnetic field was ob- 
served. The non-equilibrium conductance matches the 
previously introduced universal function of two parame- 
ters with expansion coefficients, a = 0.18 and 7 = 1.65, in 
quantitative agreement with predictions for the infinite-C/ 
Anderson model, and consistent with the allowed range 
for the finite U asymmetric Anderson model. We con- 
clude that InAs nanowires are promising new objects to 



be used in future mesoscopic transport experiments, in- 
cluding highly quantitative studies. 

There is one experimental observation, however, that 
is strikingly at odds with theoretical expectations: the 
conductance G(B) at low temperatures shows a much 
stronger magnetic field dependence than expected from 
theoretical calculations for the single-impurity Ander- 
son model (see Fig. OD^d)). As possible cause for this 
unexpected behavior, we suggest spin-orbit interactions, 
which are known to be strong in InAs nanowires^ The 
occurrence of a Kondo effect is compatible with the pres- 
ence of spin-orbit interactions, since they do not break 
time-reversal symmetry. However, they will, in general, 
modify the nature of the spin states that participate in 
the Kondo effect £2r£^ In the present geometry, where 
spin-orbit interactions are present in the nanowire (but 
not in the leads), there will be a preferred quantization 
direction (say n so ) for the doublet of local states. In gen- 
eral, n so is not collinear with the direction of the applied 
magnetic field, B. The local doublet will be degener- 
ate for B = 0, allowing a full-fledged Kondo effect to 
develop as usual in the absence of an applied magnetic 
field. However, the energy splitting of this doublet with 
increasing field will, in general, be a non-linear function 
of \B\, whose precise form depends on the relative direc- 
tions of B and n so . According to this scenario, the mag- 
netoconductance curves measured in the present work 
would not be universal, but would change if the direction 
of the applied field were varied. A detailed experimen- 
tal and theoretical investigation of such effects is beyond 
the scope of the present paper, but would be a fruitful 
subject for future studies. 
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